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Abstract 



We discuss solutions of the one-dimensional telegrapher's equation 
in the presence of boundary conditions. We revisit the case of a radi- 
ation boundary condition and obtain an alternative expression for the 
already known Green's function. Furthermore, we formulate a backre- 
action boundary condition, which has been widely used in the context 
of diffusion-controlled reversible reactions, for a one-dimensional tele- 
grapher's equation and derive the corresponding Green's function. 

1 Introduction 

Although the diffusion equation [51 3] 



can be successfully applied to a broad range of different phenomena, 
it suffers from a pathological defect: The implied signal propagation 
speed is infinite in the sense that any perturbation 5p{x, t) localized at 
one point in space instantaneously influences t) at any other point. 
Still, the diffusion equation provides an excellent approximation, as 
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long as the time scale on which inertial effects must not be neglected is 
small compared to the relevant time scale of the considered problem. 
This is not always the case, though. For instance, the movement of 
animals and biological cells proceeds with a rather well-defined and 
finite velocity [5 and the diffusion equation approximation becomes 
invalid. 

The telegrapher's equation [9l [HI [191 [13 may be considered as a 
generalization of the diffusion equation that overcomes the described 
deficiencies. Being a hyperbolic differential equation, it exhibits wave- 
like features on short time scales while diffusion-like features dominate 
on longer time scales. In particular, the signal propagation speed is 
bounded and well-defined. Like the diffusion equation, the telegra- 
pher's equation may be viewed from a macroscopic, phenomenological 
and, at least in one dimension, from a microscopic point of view. 

The phenomenological approach starts with Fick's first law that 
asserts that the probability flux is linearly related to the gradient of 
the probability density 

j{x,t) = -D . (1.2) 



Combining Eq. (1.2 1 with the conservation law 

immediately yields the diffusion equation. Consequently, Ficks' first 
law and the telegrapher's equation are not compatible with each other 



and one has to seek a modification of Eq. (1.2). A suitable general- 
ization was suggested by Cattaneo in an attempt to replace Fourier's 
law in the theory of heat waves. In the present context Cattaneo's 
equation reads as 

t'^-(o'I+A (1.4) 



dt \ dx 

where T sets the relaxation time scale. Obviously, for T — > one 



recovers Eq. (1.2). We note for later use that 



J{x,t)^^^ f e^(*-)/^^?^dr (1.5) 



satisfies Eq. (1.4), given that one identifies the diffusion constant with 
the speed c and the relaxation time T according to D = c^T. Cat- 



taneo's equation Eq. (1.4 1 combined with the conservation equation 



Eq. (1.3 1 leads to the telegrapher's equation 
d^p 1 dp 2 (^^1 



d^p , 1 dp 2 d^P /I a\ 

c (1.6) 
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Note that T interpolates between the baUistic and diffusion regime 



in the sense that Eq. (1.6 1 reduces to the wave equation in the hmit 



T — > cx), while one arrives at the diffusion equation Eq. (1.1 1 in the 
limit T ^ 0, c ^ oo, c^T ^ £) = finite. 

Turning to a microscopic point of view, we recall the well-known 
fact that the diffusion equation can be obtained as the continuum limit 
of the simple uncorrelated random walk, also referred to as Brownian 
motion. Similarly, at least in one dimension, the telegrapher's equation 
emerges from a correlated random walk that takes into account correla- 
tions in the direction of movement [71 [TSl IH] • In this context Brownian 
motion appears as a contracted description of the more general cor- 
related random walk, being valid only if the characteristic correlation 
time decays sufficiently fast. The persistence of the correlated walker, 
i.e. the tendency that each step points in the same direction as the 
previous one, necessitates to consider two separate probability density 
functions (PDF). Let a(a;,t|a:o) denote the PDF to find the particle at 
X, given that it was at xq for t — Q and that it is moving in the posi- 
tive X direction at time t, and let h{x,t\xo) the corresponding PDF for 
moving in the negative x direction. It follows that the PDF p{x,t\xQ) 
independent of the direction in which the particle is moving at t is 
given by the superposition 

p{x, t\xa) = a{x, t\xo) + b{x, t\xQ). (1.7) 

It can be shown that the time evolution of p{x, t\xo) is governed by the 



telegrapher's equation Eq. (1.6) 14 . 

The broad range of potential applications of the telegrapher's equa- 
tion [12], its blend of wave and diffusion-like features and relation to 
other important equations like the diffusion, wave and Dirac equation 
[8] , motivate a thorough study of its solutions. 



In this paper, we are interested in solutions to Eq. (1.61 in the 
presence of a single boundary [TH [TH [T3]. Previously, a radiation 
boundary condition (BC), the purely absorbing and reflecting limiting 
cases and the corresponding Green's functions (GF) have been de- 
rived jl4j . Here, we will revisit the absorbing and radiation boundary 
cases. We will obtain an alternative form of the GFs by directly in- 
verting the Laplace transform via the Bromwich contour integral and 
we will discuss how the alternative forms relate to the expressions ob- 
tained earlier. Then, we will discuss and derive a backreaction BC 
suitable for a one-dimensional telegrapher's equation. The backreac- 
tion BC generalizes the radiation BC, which only describes irreversible 
reactions between the particle and boundary: Once the particle gets 
absorbed, it stays trapped forever. In contrast, the backreaction BC 
incorporates desorption or return of population at the boundary. Con- 
sequently, the backreaction BC figures prominently in the theory of 
reversible diffusion-influenced reactions [3] and the GF of the diffusion 
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equation subject to that BC have been obtained in one, two, and three 
dimensions [21 [TTl [11] . Finally, we will focus on the purely absorption- 
desorption case and calculate the associated GF by make use of the 
Bromwich contour integral again. In his way, we can derive and discuss 
all obtained GFs in a coherent way. 



2 Radiation boundary condition 



We will start by summarizing some results obtained earlier [T?]. Eq. ( 1.6 ) 
has to be supplemented by initial conditions. A choice that can be mo- 
tivated from an underlying random walk picture is 



p{x,t = 0|a;o) 
dp 

In 



= 0. 



(2.1) 
(2.2) 



The solution to the telegrapher's equation that satisfy these initial 
conditions is known as the free-space GF [Ti] 



Pfroc(a;,i|xo) = e 



_ ,-t/(2T) 



/o(i, |x-xo|), 



where 



/o(i,2;) 



1 r/ . N Q{ct-x) 



2cT 



(2.3) 

, (2.4) 
(2.5) 



Here Q{ct — x) denotes the Heaviside step function and Iq, Ii refer to 
the modified Bessel functions of first kind [1, Ch. (9.6.)]. For later 
reference we note that the Laplace transform of fo{t, x) is 



fois,x) 



2c 



1 

4T2 



(2.6) 



The free-space GF can be used to obtain GFs that satisfy the telegra- 
pher's equation in the presence of boundaries. Typically, one seeks a 
solution of the form 



p{x, t\xQ) = Pbcc{x, t\xQ) + h{x, t), 



(2.7) 



where h{x, t) has to be a solution of the telegrapher's equation Eq. ( 1.6 1. 
Furthermore, it is chosen in such a way that Eq. (2.7) satisfy the spec- 
ified BC. Henceforth, we will assume without loss of generality that 
the boundary is located at a; = and that a; > 0, a:o > 0. 
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In Ref. [H] it has been shown that the appropriate radiation BC 
is quite different to the one in the diffusion case and that it takes the 
form 



dp 

dx 



x=0 



1 

T 



x=0 



Here, /3 denotes the probabihty that a particle reaching the boundary 
gets actually absorbed 



a{x ^ 0, i|xo) = (1 - I3)b{x 0, t\x^). 



(2.9) 



Thus, /3 = 1 corresponds to the purely absorbing limit, while /3 = 
corresponds to a purely reflecting boundary. Furthermore, seeking a 
solution of the form 



p(x,t|xo) = e-*/(2^)p(x,i|a;o), 
one finds for the Laplace transform of P{x,t\xi^) [T3] 

g -p\x-xa\/c _ ^^-p{x+xo)/c 



P{x,s\xq) 



1 r 

2T 



2cp 



where 



P 



4T2' 



2T ' 



(2-/?)p 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



It is easy to show that for purely absorbing and reflecting BCs the 
expression for Vl reduces to 



n 



2T(s-p), 
1, 



if/? = l, 
if /3 = 0, 



(2.14) 
(2.15) 



respectively. Taking into account these identities, Eqs. (2.6), (2.11) 
and elementary properties of the Laplace transform, one arrives at 
the following form for the GFs corresponding to purely absorbing and 
reflecting BC 



p,hs{x,t\xn) - e-*/^^[fo{t,\x - xo\) - hit,x + xn)], (2.16) 
Prefix, t\xa) = e-'^^^[fo{t,\x-xo\) + fo{t,x + xo)], (2.17) 

where the function fi is given by 
e{ct - x) 



8cT 

ct — X 
ct + X 



/o(w) + 2 

I2{U) 



ct — X 
ct + X 



1/2 



Iiiu) + 



(2.18) 
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Finally, to obtain the inverse Laplace transform in the case of a radi- 
ation BC one can proceed by expanding the denominator in a power 
series to find the following expression for the GF satisfying the BC 

Eq. p:8| m 



•i,hs{x, t\xo) + (1 - /3)prcf(2^, t\xo) + 

^ OO 

— /3(1 - /3)e-*/2^ J2 C"5n(i, X + xo), (2.19) 



8cT 



n=0 



The coefficients are defined by 

Co = 1, 

Cl = 4-/3, 

C2 = (2-/3)(3-/3) + l, 

cn = (2 - /3)3(1 - /3)"-^ n>3, 



gnit,x) = 0{ct-x) 



ct ~ X 
ct + X 



n/2 



ln{u). 



(2.20) 



2.1 Alternative expression for the Green's function 

In what follows, we will take a different route to derive alternative 
expressions for the GFs subject to purely absorbing and radiation BC. 
By making use of the Bromwich contour integral, we will derive an 
alternative form of the GFs as an important preparatory step. In fact, 
by employing for all discussed BC the same technique, hardly any 
extra work will be necessary to establish the expression for the GF 
that satisfies the backreaction BC. 

Focusing first on the purely absorbing case, we start again from the 
Laplace domain and write 



Pabs(a::, s|a;o) = Pbcc{x, s\xq) + Khs{x, s). 



(2.21) 



Now, our strategy will be to invert /labs directly by using the Bromwich 
contour integral. Henceforth, to streamline the notation, we define 



X + Xq 



c 



Then, taking into account Eqs. ( 2.11 1, ( 2.12 1 and ( 2.13 1, the Bromwich 
contour integral takes the form 



{x,t\xn) 



1 

27ri 



7+ioo 



(x, s)ds 



7 + 200 



7— zoo 



ds. (2.22) 



6 



To proceed we have to discuss the analytical structure of the integrand, 
see Fig. [T] First, we note that it has two branch points, at s = — 1/T 
and s = 0, necessitating the specification of two branch cuts. We 
choose both branch cuts along the negative real axis, leading to a 
potentially discontinuous integrand in the region s < 0. Therefore, we 
employ the keyhole contour depicted in Fig. [l] It turns out that the 
branch cuts' effects cancel each other along the overlapping region s < 
— 1/T, i.e. along C_^'^, CZ'^jrp, and the integrand is only discontinuous 
for — 1/T < s < 0. Therefore, the sole contributions coming from the 
integration on each sides of the branch cut are due to the integrals 



along C_^^y and Cq 



1/T 



Along C^i^j. we choose s = re*^ and it follows 
s + a = \J a — r and \/s = i\pr. We get 



C'li/^ Jo i^Jr{a-r) ^a-r + i^r 

(2.23) 

— 1/T 

Similarly, to evaluate the integral along Cq , we choose s = re~^'^ 
and note that the ensuing integral is the negative complex conjugate 



of the integral given in Eq. (2.23) 



e^*/iabs(a;, s)ds 



(2.24) 



where * denotes complex conjugation. Furthermore, the integrand has 
no poles neither inside nor on the contour and the contributions from 
the small circles Ce2,Cej^,Ce^ around the origin and —1/T, respectively, 
vanish. Thus, using Eqs. (|2.24[), (|2.23[), (|2.22| and making the substi- 



tutions r — >■ </? + l/2a, ip — > l/2a^, we finally arrive at 



aB(ci-fc)^_i/2„t 



47rc 

cos{l/2kay/l-e 



^-l/2atC 



(1-0 



U + sm{l/2kax/l-e) 



dC (2.25) 



The question arises how the found expression relates to the expres- 



sion obtained in Ref. [Tlj, cp. also Eqs. (2.16) and (2.18). Actually, 
one can explicitly show that Eq. (2.25) can be cast in an integral free 



form. To this end we first note the identity 

1 ,-i/2a., cos(l/2afcyr 3g)^^ ^ 4(l/2ay^^). (2.26) 

TT J -I ^\ - ^2 



Then, it follows that the terms involving the cosine in Eq. (2.25) can 



be rewritten by taking the first and second order derivative of the 
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expression on the rhs of Eq. (2.26) with respect to t. Also, upon 



integration by parts we can bring the term involving the sinus in a 
form that can be expressed analogously. 

The general case of a radiation boundary conditions can now be 
treated in a similar way. Again, the part of the GF that takes into 
account the BC can be expressed by a Bromwich integral 



/irad(a;,t) 



1 



Attic 



7+ioo 



s + a (1 — ri)^ s + a — (1 + iii)^/s 
sjs{s + a) (1 - -nWs + a + (1 + r^)^s 



ds. 



(2.27) 



where rj = \ ~ 13. Because the analytical properties of the integrand 



in Eq. (2.27) are analogous to the case of purely absorbing BC, we 



employ the same keyhole contour and proceed in the same way to find 



47rc 



1 1 + 772 + 2?7^ 



(i-e)x 



^°^^'^^[2. + (1 + ,^),] + [1 - sin(^yr^) 



d^. (2.28) 



Here, two remarks are in order. First, as 77 — > 0, one recovers the 
expression obtained for /labs Eq. (2.25). Second, we may expand the 
denumerator 1 + 77^ + 2ry^ in powers of 2ri^/ (1 + 77^), because 277/ (1 + 
77^) < 1 for 77 < 1. As explained for /labs, the resulting integrals can be 
expressed as suitable time derivatives of the expression on the rhs of 



Eq. (2.26), and in this way one may recover the series expansion given 



in Eqs. (2.19), (2.20) 



3 Backreaction boundary condition 

Next, we will turn to reversible interactions with the boundary: when 
the particle gets trapped this does not mean necessarily that it stays 
trapped forever. First, we will derive a backreaction BC that is suit- 
able for the one-dimensional telegrapher's equation. The derivation 
will follow the line of reasoning in Ref. 14J and it will turn out that 
the backreaction BC is somewhat more complicated than the corre- 
sponding BC in the case of diffusion, resembling the situation for the 
radiation BC. 

Conventionally, the salient assumption underlying the backreaction 
BC is that the rate of desorption is proportional to the total absorbed 



population. Hence, generalizing Eq. (2.9), wc postulate that at the 



boundary the PDFs that take into account the directions are related 
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by 



aix,t\xo)U=o = (1 - l3)bix,t\xo)U^o + " S{t\xo)]. (3.1) 
Here, S{t\xQ) denotes the survival probability, defined by 

nOO 

S{t\xo) ^ / p{x,t\xo)dx, (3.2) 



and therefore 1 — S{t\xo) is the total absorbed population. The total 
probability density becomes 



p{0,t\xo) = a{0,t\xo) + b{0,t\xo) 

= {2~mO,t\xo) + ^[l~S{t\xo)]. 



Additionally, we can conclude that 



a{0,t\xo)-biO,t\xo) = - 



2-/3 



p{0,t\xo)- k[1- S{t\xo)] 



Moreover, one has the dynamic equation |14j 



dp 
dx 



d 1 



(3.3) 



(3.4) 



(3.5) 



which is valid in particular at the boundary a; = 0, and hence we arrive 
at the backreaction BC for the one-dimensional telegrapher's equation 



dp 
dx 



13 (dp 



2- p\dt 



dS{t\xo) 

' dt 



1 

T 



p\x=o- S{t\xo)) 



(3.6) 

To shed more light on the nature of the parameter k it is instructive to 
consider the limiting case of the BC Eq. ( 3.6 ). Taking the limit c — >■ oo, 
T ^ and ^/(2 - ^) ^- such that c^T =: D and /3/(2 - P)c =: Ka 
are kept finite, we recover the radiation-backreaction BC ^ for the 
diffusion equation 



D 



dp 

dx 



KaP\x=0 - Kd{l ~ S{t\xo)), 



(3.7) 



£C = 



where := KqK. On the other hand, when we take the same limit 
as before with the exception that f3/{2 — f3)c — ?> oo, we arrive at the 
absorption-backreaction BC [2] 

p\x=o = 4l-S{t\xo)]. (3.8) 

To derive a GF of the telegrapher's equation subject to the BC 
specified by Eq. (3.61, one has to find a suitable expression for the 
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survival probability. To this end, we briefly recall how one proceeds 
in the case of diffusion. From the equation of continuity Eq. (1.31 and 
the definition of the survival probability Eq. (3.2 1 it follows 

dS{t\xo) 



dt 



: j{x,t\Xo)\x=0. 



(3.9) 



Assuming that the current vanishes at infinity and that Fick's first law 



Eq. (1.2 1 is valid one obtains 
S{t\xo) = 1 - 



dx 



dt'. 



(3.10) 



x=0 



Obviously, Eq. ( 3.10 1 cannot be used in the context of the telegrapher's 
equation, because we already discussed that Fick's law has to be re- 
placed by Cattaneo's equation Eq. (1.4). Hence, we combine Eq. (1.5) 
and Eq. (3.9), leading to 



S{t\xo) = 1 



"'0 



OX 



(3.11) 



The derived expression for the survival probability can be intro- 
duced in the BC Eq. (3.6 1. Henceforth, we will focus on the case 

(3.12) 



(3 = 1. Then, in the Laplace domain the BC becomes 



Kc^ ^ dp 
dx 



x=0 



As before, we make the ansatz 

p{x, S\xa) = Phcc{x, s\xq) + ft.back(2^, s) . 



(3.13) 



The function ft.back(a^, t) has to be a solution of the telegrapher's equa- 
tion and hence takes in the Laplace domain the form 



(3.14) 



The constant A{s) must be chosen in such a way that Eq. ( 3.13 1 satisfies 
the BC. We obtain 



A{s) 



1 s + a g + g- (l + f)Vg(g + a) ^_5c,/;(^ 
2c ^s{s + a) s + a + (1 + ^) ^ s{s + a) 

(3.15) 

As before, we can calculate /^back(3:,t) via the Broniwich integral 



''back 



(x,<) 



1 



S + a 



^/s{s + a) 



gts~k^s(s+a) ^ 

7 — iCJO 

y^T^- (1 + f )ys 



+ a + (1 + ^ 



ds. 



(3.16) 
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We can closely follow the line of reasoning presented for the case of 
absorbing and radiation BC. We may specify the same branch cut 
structure, employ the same integration contour and make the same 
substitutions. Thus, we arrive at 



, , , aQ(ct — k) 1 



Aire 



1 c^K^ + a/2{a ~ 2ck){1 + 



cos{l/2ka^/T^] 



Here, we introduced the functions 



n(0 + sm{l/2ka^l-e)r{0 



m 



1 



1 



,2„2 



a ^ + -a(a — 2ck)^ — ck + c k 
-a^^ + 2'^('^ ^ 2ck). 



(3.18) 
(3.19) 



In the limit k — > 0, we recover the result for the purely absorbing BC 



Eq. (2.251. Furthermore, we note that 

f (a - 2ck) 



1{a-2cK) + c^K^ 



< 1, 



(3.20) 



and therefore in this case also it is possible to expand the denumerator 
in a power expansion in 



f (a - 2ck) 
I (a - 2ck) + c^K^ 



Using Eq. (2.26) this means that we can dispense with the integral 



representation altogether in the case of the backreaction BC also. 



(3.17) 
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